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While biogeochemical and physical processes in the Southern Ocean are thought to be central to 
atmospheric CO2 rise during the last deglaciation, the role of the equatorial Pacific, where the largest CO2 
source exists at present, remains largely unconstrained. Here we present seawater pH and PCO2 variations 
from fossil Pontes corals in the mid equatorial Pacific offshore Tahiti based on a newly cahbrated boron 
isotope paleo-pH proxy. Our new data, together with recahbrated existing data, indicate that a significant 
pCOa increase (pH decrease), accompanied by anomalously large marine "C reservoir ages, occurred 
following not only the Younger Dryas, but also Heinrich Stadial 1 . These findings indicate an expanded zone 
of equatorial upwelling and resultant CO2 emission, which may be derived from higher subsurface dissolved 
inorganic carbon concentration. 

Understanding the past condition of the surface ocean carbonate system and air-sea CO2 exchange is crucial 
to projecting future changes in the carbon cycle under ongoing anthropogenic global warming and ocean 
acidification. Atmospheric CO2 concentration increased by as much as 80 \iatm during the last deglacia- 
tion, with —50 |iatm released during Heinrich Stadial 1 (HSl, from 17.5 to 14.6 kyr), followed by an additional 
—30 |j,atm during the Younger Dryas (YD, from 12.9 to 11.7 kyr) (ref 1). While the Southern Ocean is generally 
considered to be central to the deglacial CO2 rise^ the contribution from other oceanic regions remains 
relatively uninvestigated'" Information on the partial pressure of CO2 (PCO2) is needed to directly constrain 
past air-sea CO2 exchange, and this can be reconstructed from boron isotopes (8"B), a marine carbonate pH 
proxy'''''^ Regions where surface seawater CO2 is out of equilibrium with the atmosphere are ideal for such 
studies, and the equatorial Pacific is particularly well suited because it represents the largest global CO2 source at 
present (e.g. ref 16). Corals broadly distributed in tropical to subtropical areas constitute excellent high-resolu- 
tion geochemical archives for paleo-C02 studies because they may be precisely radiogenically dated (U-series)^ '', 
unlike foraminifers that are affected by uncertainty in radiocarbon reservoir age (R). 

Two previous studies have attempted to constrain equatorial Pacific CO2 changes. Palmer and Pearson'" 
showed increased CO2 emission during the last deglaciation in the western equatorial Pacific (WEP) from 
5"B measurements on the planktonic foraminifer {Globigerinoides sacculifer) in a sediment core recovered from 
offshore Papua New Guinea (ERDC-92, Fig. la). Further east, DouvUIe et al." performed S"B on fossil corals 
from the Marquesas (9.5°S 139.4°W, Fig. la) and also demonstrated increased CO2 release at the end of the YD. 
However Douville etal." did not observe a significant CO2 release during HSl, complicating interpretation of the 
equatorial contribution to deglacial atmospheric CO2 rise. Integrated Ocean Drilling Program Expedition 310 
(lODP Exp. 310)'" drilled the outer reef slope at Tahiti (17.6°S 149. 5°W, Fig. la) recovering fossil corals from an 
open ocean environment spanning HSl, which enable us to assess the issue. 

Results 

8"B-pH calibration. This study establishes a new empirical S"B-pH calibration utilizing, for the first time, 
anthropogenic ocean acidification. Empirical calibration is needed to overcome the observed offsets from a 
theoretical 5"B-pH curve in culture experiments for zooxantheUate corals'* '^. (pH is reported using the total 
hydrogen scale, hereafter pH for simplicity). There are two primary approaches to overcoming the reported 
offsets. One is an empirical approach''* that assumes constant offsets in measured and theoretical §"B ("offset a;" 
see Methods), while the other is an observational approach'^ that considers potential pH modification by 
calcifiers. 
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Figure 1 | Sea surface pC02 in the equatorial Pacific and vertical sections of the South Pacific, (a) Equatorial Pacific locations discussed in the text with 
sea surface pCOj for the reference year AD 2000 in which red and blue represent CO2 sources and sinks, respectively'". Vertical sections from 60 'S to 30 'N 
of (b) oxygen, (c) silicate, (d) DIG, and (e) A^C. All data are downloaded and plotted using Ocean Data View software, version 4.5.5 (ref 60). Inverted 
black triangles show locations of Tahiti and Marquesas. The white star and diamond in e indicate the locations of cores MV99-MC99/GC31/PC08 
(refs. 36,43) and V21-31 (ref 42), respectively. 



The potential for pH-modification is of great concern for S"B- 
based reconstruction of pH due to the implications for atmospheric 
PCO2 calculation. Such a phenomenon is consistent with indirect pH 
measurements of internal calcification fluid using pH sensitive dye 
that suggests a higher pH than ambient seawater, creating better 
conditions for calcification". A 'ApH' concept in the 5"B-pH cal- 
ibration that reflects pH differences in seawater (pHsw) and internal 
calcification fluid has been proposed'^, however usage of this pro- 
posed relationship to calibrate 8"B of modern Porites spp. from 
Tahiti and Marquesas resulted in unrealistically high values (e.g. 
—8.34 in AD 1991), well above reported estimates (e.g. refs. 20,21). 
Therefore the present study employs the empirical equation" here 
(see Methods). 

Ocean acidification was estimated from a combination of in situ 
fC02 values in the surface ocean, atmospheric CO2 concentration 
directly measured at the Mauna Loa observatory in Hawaii since AD 



1960 (ref. 22), and CO2 concentrations within bubbles trapped in an 
Antarctic ice core^' (Figs. 2, S1-S4; Supplementary Methods). These 
data were then fit to the previously reported S"B measurements of 
Porites spp. (refs. 11,24), which are for the years AD 1991, 1950, and 
1700 (Fig. 2, see Methods and Supplementary Methods for details of 
the 5"B-pH calibration and pH estimation since the Industrial 
Revolution). 

pH and PCO2 reconstruction. Using our revised calibration, we 
reconstructed pH from our new 5"B measurements on Tahitian 
corals, as well as from previously reported data" from both the 
Marquesas and Tahiti, and the overall result is consistent with the 
WEP foraminifer 5"B variations'" (Fig 3a and b). The oldest coral 
sample, dated to 20.7 ka BP during the last glacial maximum (LGM), 
exhibits a relatively high pH (8.26). From 15.5 to 9.0 ka BP, pH is 
generally constant within uncertainty (8.15-8.22) and consistent 
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Figure 2 | Measured pC02 from Mauna Loa^^ and Law Dome ice core^' 
(red and blue lines) and estimated pH variation at Tahiti (black line). 

Green triangles are newly calibrated pH from previously reported 5"B of 
modern corals from Tahiti and Marquesas"'^* (Table SI). Estimated pH 
from the AD 1700 Marquesas coral was scaled by 0.04 to correct for offset 
from Tahitian coral values. Error bars are 2a. 

with the preindustrial value of 8.20. Four notable pH excursions are 
associated with HSl and the YD. Two of our samples exhibit 
anomalously low pH at the end of HSl (8.13 at 15.15 ka and 8.09 
at 14.99 ka BP), in addition to those at end of the YD at the 
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Figure 3 | Deglacial 8"B, pH and pC02 variability in the equatorial 
Pacific, (a) Reported S"B values of planktonic foraminifera with 2o 
uncertainty from ERDC-92 (ref 10). Age model is from the original 
publication, (b) In situ pH reconstructed from 8"B of fossil Pontes spp. 
using our new calibration. Red and green circles are from Tahiti by this 
study and Douville et aV\ respectively. Blue circles are from Marquesas by 
Douville et after correction by +0.04 pH units, (c) Derived pC02 of 
surface water around equatorial South Pacific Ocean (same colors as b) 
and atmospheric pC02 on the GICC05 timescale'. All error bars are 2a. 
YD, Younger Dryas; B/A, BoUing/AUerod; ACR, Antarctic Cold Reversal; 
HSl, Heinrich Stadial 1; LGM, Last Glacial Maximum. 



Marquesas". The low pH following HSl had been previously 
undetected at this location. Calculation of pC02 (see Methods) 
reveals deglacial values significantly above those of the atmosphere 
(Figs. 3c and 4a). Conversely, ApC02 during last glacial and the early 
Holocene was nearly zero, suggesting air-sea CO2 equilibrium. 

Results from a different portion of the same 14.99 ka BP coral 
sample deviate by as much as 1.4%o, which corresponds to 0.11 in 
pH and 100 \iatm in pCOj (310-M0024A-llR-lW_77-90 and 310- 
M0024A-llR-lW_60-75, Table SI). Considering the average ~4 
year temporal resolution of each sample, these excursions occurred 
abruptly and persisted for several years, which differs from modern 
observations that show no clear interannual or decadal variability 
(Supplementary Methods, Fig. SI). This enhanced variability, which 
is also observed in Sr/Ca derived SST results from another Porites 
colony recovered from lODP Exp. 310 (ref. 25), may relate to Tahiti's 
location at the rim of equatorial upweUing cell (Fig. la). Taken 
together, PCO2 (pH) records indicate that the equatorial Pacific 
became a larger CO2 source during the last deglaciation with excur- 
sions at the end of HSl and the YD. 

Marine "C reservoir age compilation. Compiled marine "C 
reservoir ages (R) throughout the equatorial Pacific resemble PCO2 
variability (see Methods; Figs. 4 and S5; Table S2). Larger and more 
variable values of R are evident in Tahiti during HSl and the YD, and 
enhanced R variability is also seen in the Marquesas (Figs. 4 and S5). 
Paterne et al?'' sub-sampled different parts in the same fossil coral 
skeleton and analyzed both "C and U/Th. They observed no 
difference in U/Th dates, but a much larger difference in "C. 
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Figure 4 | Differences between modern and deglacial PCO2 and R 
(ApC02 and Rdiff). The last deglaciation is characterized by larger and 
more variable PCO2 and R values, (a) Atmosphere-ocean ApC02 based on 
Tahiti and Marquesas corals and Antarctic ice core'. Figure legend is same 
as Fig. 3. (b) Rdin around the equatorial Pacific Ocean calculated from fossil 
corals (red: offshore Tahiti; orange: reef crest of Tahiti barrier reef*, blue: 
Marquesas^'; light blue: Kiritimati"^'; black: Mururoa^") . Horizontal dashed 
lines represent ApC02 = 0 and R^itj = 0. 
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Possibilities of either a diagenetic alteration or a change in R were 
suggested. The latter is more probable because a large variation in R is 
also suggested from Vanuatu coral at 11.7-12.4 ka (—400 years; 
during the YD, ref 27). 

Though reported R^ig (difference between calculated R and mod- 
ern R; see Methods) data around the upweUing zone during the LGM 
are sparse, calculations with the new Lake Suigetsu datasets^** sug- 
gests no substantial change in R (see Methods; Figs. 4 and S5). This 
implies that the CO2 exchange rate in the surface equatorial Pacific 
during the last glacial was almost the same as present, which supports 
the above-mentioned observation that ApC02 is essentially equival- 
ent to zero and indicates that anomalous R values are limited to the 
last deglaciation. 

Discussion 

PCO2 variability in subtropical oligotrophic water can be explained 
by mixing of water masses that exhibit distinctly different dissolved 
inorganic carbon (DIG) concentrations. A southward migration of 
the intertropical convergence zone (ITCZ) that partly controls ther- 
mocline depth is hypothesized during Heinrich Events including 
HSl and the YD (e.g. ref 29). At present, the ITCZ does not seem 
to affect surface pC02 variability (Fig. la), and if it is displaced 
southward, the locus of equatorial upweUing remains at the equator 
due to the influence of inter-hemispheric asymmetry of Coriolis 
force (e.g. refs. 30,31). Enhanced upweUing (shallower thermocline. 
La Niiia-like conditions) or increased subsurface DIG concentration 
are more likely to drive PGO2 variability based on sedimentary evid- 
ence from the equatorial Pacific for higher nutrient content, e.g., 
enhanced biogenic opal export production and lower stable carbon 
isotopes (5''G) (TT013-PG72, ODP Site 1240 and TR163-19)'"= 
(Fig. la). Semi-conservative radiogenic neodymium isotopes (sNd) 
from sediment cores at the eastern equatorial Pacific (EEP) (ODP 
Site 1240) and off Baja Galifornia (MV99-MG99-GG31/PG08) indi- 
cate stronger subsurface water transport from the south'''''"' (Fig. 1). 
Govariation of geochemical properties between the Southern Ocean 
and the equatorial Pacific suggest a subsurface connection during the 
last deglaciation (e.g. refs. 32-38). Thus, PGO2 variability may be 
explained by an increase in DIG in the upwelled, subsurface water 
masses as opposed to physical processes alone. 

Water mass subduction along the subantarctic front, mainly off 
GhUe^'^, forms Subantarctic Mode Water (SAMW) and Antarctic 
Intermediate Water (AAIW) that upwells at the equatorial Pacific 
via the Equatorial Undercurrent (EUG) (Fig. 1). SAMW and AAIW 
are characterized by higher/lower concentrations of oxygen/sUicic 
acid (Fig. lb,c). It is suggested that the abyssal DIG reservoir around 
the Southern Ocean increased during the last glacial period"" 
which would have contributed to lower atmospheric pG02. 
Garbon dioxide was released to the surface through deep ocean vent- 
ilation during HSl and the YD (refs. 2-9), however export produc- 
tion was insufficient to fully compensate the increased carbon flux*". 
This is consistent with residual radiocarbon content (A'^G) of inter- 
mediate water at the EEP (V21-30)''' and off Baja Galifornia (MV99- 
MG99-GG31/PG08)'" that indicates anomalously older water was 
incorporated into SAMW/AAIW (Fig. 1), as well as with depleted 
8"C of surface and lower thermocline dwelling foraminifers from 
sediment cores at both equatorial (TT013-PG72, ODP Site 1240 and 
TR163-19) and South Pacific sites'^ '^'''''"'. Moreover, enhanced 
export production of biogenic opal suggest more silicic acid was 
transported via the EUG to thermocline water at the equatorial 
Pacific (V19-30 and TT013-PG72) without being consumed com- 
pletely within the Southern Ocean'*'''''". Stronger Ekman transport in 
association with sea ice retreat and a poleward shift of southern 
westeries is suggested to be a driver''''. 

A similarity between R and PGO2 variability during the last degla- 
ciation supports an interpretation that older DIG was incorporated 
to subtropical surface water through mixing with SAMW/AAIW, 



though, contrary evidence comes from the current formation sites 
off GhUe*'' and New Zealand*. However, a key sediment record off 
Ghile was recently reevaluated, and the new interpretations indicate 
stronger upweUing and subsequent larger R in surface water in the 
Southern Ocean', which agrees well with our interpretation. Yet, 
further work is stiU needed to fully understand both the physical 
and biogeochemical dynamics in the Southern Ocean and the equat- 
orial Pacific"'. 

Positive ApG02 indicates GO2 flux from the ocean to the atmo- 
sphere. Previous studies'"'" indicated that the equatorial Pacific con- 
tributed to deglacial GO2 rise, however the timing of anomalously 
higher PGO2 events recorded in radiogenically dated fossil corals do 
not systematically correspond to those of atmospheric GO2 rise 
recovered from Antarctic ice core on the GIGG05 timescale' 
(Figs. 3 and 4). Moreover our new calibration reveals a modest 
GO2 emission continued through the Boiling/ AUerod/ Antarctic 
Gold Reversal when no atmospheric GO2 increase is observed 
(Figs. 3 and 4). Though we demonstrate that the equatorial Pacific 
became a larger GO2 source during the last deglaciation, it is too early 
to conclude its exact contribution to atmospheric GO2 rise. The 
Southern Ocean is suggested to be central in GO2 degassing'' ' and 
the contribution of the terrestrial biosphere should be further eval- 
uated"". More evidence spanning the YD and the early part of HSl, in 
particular the sharp rise in atmospheric GO2 and the sudden drop of 
8"G of GO2 (refs. 1,46), as weU as more spatial coverage is needed. 

Methods 

6"B analyses. The S^^B values of fossil Pontes spp. were measured following the 
protocol of Ishikawa and Nagaishi"*^. Fossil corals were screened for diagenetic 
alteration with X-ray diffraction and geochemical analyses, as well as visual using a 
Scanning Electron JMicroscope'^'*. Bulk sampling was conducted along the growth axis, 
and time resolution of each sample is several years {1-8 years) depending on growth 
rate of each coral"^**. Typically 6 mg of carbonate was used for 5^^B measurement. After 
removals of organic matter using 30% H2O2 for —12 hours, boron was purified by 
cation and anion exchange using AG 50 W X12 and 1-X4 resin (Bio-Rad, USA) and 
then 5^'B were measured using the positive polarity thermal ionization mass 
spectrometer (P-TIMS; Thermo Finnigan TRITON) installed at Kochi Core Center, 
lapan Agency for Marine-Earth Science and Technology. All reported 6^'B values are 
the mean of duplicate analyses (Table SI). Repeated analysis of the JCp-1, carbonate 
standardprovidedbyGeologicalSurvey of Japan yielded 24.21 ± 0.18%o (2cj, n = 18), 
which is the finest precision to date"*^. Differences between the duplicates are 0.08%o 
on average with the largest one of 0.28%o {Table SI), which is within the measurement 
uncertainty of JCp-1. We conservatively report ±0.18%o as the analytical uncertainty 
of our 6^'B measurements. 

8^^B-pH calibration and PCO2 calculation. First, the analytical procedure-specific 
isotopic offset''^ was corrected using the equation bellow {modified after Zeebe & 
Wolf-Gladrow^°) in order to reduce inter-laboratory offsets in reported 5^^B values, 

clln _ sIIr , 'S"Bsw-aacked + l-0' 

O ^^carbonate-corrected — <^ ^SW-stacked "t" ~7u~Z TTT , , 

d Bs„_„casured + 10' (l) 

* ("^ S carbon ate— measured ~~ ^ -^SW— measured ) 

where 5^'Bcarbonate-corrected IS the boron isotopic value of carbonate after correction 
and S^^Bsw-stacked is the global average 5^^B of seawater. We used the recommended 

value of 39.61%o for S'^BsW-stacked (ref. 49). 6 "B^^rbonate- measured and 8"BsW-measured 

are 5'^B of carbonate samples and of seawater measured in different laboratories. 
Without performing inter-laboratory correction of JCp-1, results of Douville et al}^ 
and this study differ by 0.25%o, equivalent to —0.03 pH units. After the correction 
using the reported value of 5^'Bsw-measured (refs. 11,47), the difference improved to be 
0.07%o, indicating that correction should be performed. 

After the correction, the empirical calibration equation reported by Honisch et al}* 
was utilized. In this equation, the vital offset 'a' and fractionation factor 'a3_4' that 
yield the lowest e^ms (root-mean square error) is determined. 

-tj 1 f ~(^^^Bsw-stacked "^^^Bcarbonate-corrected "^3) \/',\ 

pH = pKB— logi (2) 

\(5 By w - stacked — ^;3- 4 * Bcarbonate-corrected + a) — 10^ * (^3-4 — 1)/ 

Three previously reported 5^^B values for modern Pontes spp. from Tahiti and 
Marquesas were fit to estimated pH at AD 1991, 1950 and 1700 (refs. 11,24) (Fig. 2, 
Table Si). In the calibration 'a' was determined as ' — 6.0' when we chose recently 
reevaluated a3_4 of '1.0272' (ref. 51). When this calibration was conducted, both 'a' 
and 'o(3_4' were allowed to vary freely, and the resultant a3_4 with the lowest e^^s was 
very close to that of Klochko et al.^^ rather than the previously accepted value of 
Kakihana et aV^ Thus we used Klochko's fractionation factor. pH was calculated from 
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both previously reported and newly obtained 6^^B values considering a local differ- 
ence of seawater pH between Tahiti and Marquesas (see Supplementary Methods for 
details). 

PCO2 was further calculated from obtained pH values from 5^^Bcarbonate-corrected 
using C02SYS program under assumptions of constant temperature, salinity and 
total alkalinity {see also Supplementary Methods for details). The values of G. 
sacdifer^° are not included here due to large uncertainties in the 5^^B-pH calibration 
(e.g. refs. 12,15). 

R compilation. Published ^^C (radiocarbon years) and U/Th ages of fossil coral 
samples obtained during lODP Exp. 310 were compiled in order to calculate residual 
radiocarbon activities (A^^C) and R. We verified via lODP sample ID and core 
photographs^^ that the exact same samples were selected (Table S2). In some cases 
different portion of the skeleton of the same coral was dated. Given that lifetimes of 
coral are generally less than several decades, temporal gaps derived from sub- 
sampling are negligible in calculations of A'*C and R. We did not use ages from 
either microbialite (carbonate created by bacteria) or encrusting coralline algae from 
equivalent down-core depths due to a possibility of post-depositional growth (for 
details, see ref 53). Calculation was done according to equations (3) and (4) where 
^*C-age is an original radiocarbon data before a local R correction^*'^^. 

I e 8033 \ 

A"C„„i,„=^^-l *10' [%.] (3) 

\ e 8266 / 



R — 8033*lnl -75- — Cage(marine)- Cage(atm) (4) 

Atmospheric A^'^C (A'^C^un) from INTCAL09 (ref 56) was used to calculate R. 
Samples that span 29-30 ka BP were calculated using the recently published Lake 
Suigetsu A^^Cg^jn dataset^^, since as much as "-100%o offsets are observed between 
INTCAL09 in this interval (Fig. S5a). Southon et aU^ points out that the INTCAL09 
curve heavily relies on Cariaco Basin varve sediment data beyond the dendro- 
chronological limit (>12.55 ka BP) and has deficits during the YD and HSl due to 
changes in local R in the equatorial Atlantic. Specifically, a difference during 15.5- 
17.0 ka BP is apparent, however it does not affect our calculation since there are no 
samples that span this interval (Fig. S5). Other and U/Th datasets for fossil corals 
from the equatorial Pacific islands were also used to calculate R including Tahiti^^, 
Marquesas^^ Kiritimati^^ and Mururoa^^ (Figs. 1,4, and S5b). R^iff denotes differences 
between calculated R and modern R (235 + 110 for Tahiti'^ 390 ± 60 for 
Marquesas^^ 335 ± 100 for Kiritimati^^ 300 ± 100 for Mururoa^^), thus it differs 
from AR that conventionally represents local reservoir age. We estimate that the 
accumulated uncertainty in the Rjifl- calculation are the sum of errors in ^^C dating, U/ 
Th dating, and modern R (Fig. S5b; Table S2). 
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